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Previous work included solving the small perturbation equation for 
nonlifting steady flow, in parameter space, using the Alternating Direction 
Implicit method (ADI).^ These solutions were obtained on a coarse grid 
with 10 points on the airfoil (Ax = . 1 ). This relatively large grid 
spacing in the region near the airfoil leads to some irregularities in 
the computed surface pressure coefficients (Cp) near the airfoil leading 
and trailing edges. To alleviate this problem, it was necessary to reduce 
l^e grid spacing, in the vicinity of the airfoil. 

In order to obtain a better definition of the flow in the region 
near the airfoil, a coordinate stretching was introduced. The stretching, 
which was Introduced by Carlson , is symmetrical and divides the flow 
field into three regions. Those regions are illustrated in Figure 1 . 

The relationships between the stretched and unstretched coordinate systems 


are 

z = a.jtan(|- n) (1) 

Xj = -x^ + agtanCjC? + 54)] + (2a) 

Xii = ^(^4 + (2b) 

Xiii ** + a2tan[|(5 - ^)] + a3tan[|(c - ^4)^] (2c) 


where a.|, a^, a^, a^, b, X4, and are constants which determine the 
geometry of the stretched coordinate system. The effect of the coordi- 
nate stretching is to increase the number of grid points near the airfoil 
and to map the infinite physical plane, -co £ z -» < x £ onto the 
finite computational plane, -1 £ n 5. 1 > -(1 ■*" C4) 5 C £ (1 ■’‘54)* 

The steady-state small perturbation equation has been approximated 
by finite differences, in stretched coordinates, and solved using column 
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relaxation. To help in choosing a relaxation factor (u), the maximum error 
(E) in the entire flow field, after 250 iterations, has been plotted as 
a function of w. E is defined as 


E = 



0 ) 


( 3 ) 


k "fch 

where „ is the value of g at the k interation. Figure 2 summarizes 
^n ,m 

E versus oj. 

Results were obtained for subcritical flow about a nonlifting para- 
bolic arc airfoil at free stream Mach numbers (M^) of ,806, .825, .85, 

k+1 k 

and .9. The solutions were considered converged when 9„ m - 9„ m was less 

n a ^ n 

-5 

than 2*10 everywhere in the flow field. The results are summarized in 
Table I. Figure 3 shov/s computed values of C at = .825 with the 

r 

airfoil thickness ratio (t) ranging from .01 to .06. 

Current efforts include the calculation of steady flows with shock 
waves and the calculation of two-dimensional unsteady flows. 

Using the method of parametric differentiation and the integral 
equation method, a formulation of the small perturbation, unsteady, non- 
linear transonic flow everythin, shockless airfoils has been developed. 
This formulation allows for a calculation of the flow field variables 
as well as the surface pressure distribution. Singularities of the 
kernel functions of both the steady and unsteady solutions have been 
removed at all points in the flow field save the leading and trailing 
edges. The integral equations for the steady and unsteady components 
of the solution for g(x,t) have been written in such a manner that a 
classical iteration procedure (Picard’s method) at each t - level is 
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appropriate. This means that the typical unsatisfactory approximations 
of the variables constituting the integrands may be avoided. 

Additional analysis to complete the evaluation of remaining integrals 
are in progress. 
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TABLE I 


SUMMARY OF RESULTS 
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Figure 1. Stretched Coordinate System 






Figure 2. Maximum Error versus Relaxation Factor 
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